APPLICATION FOR UNITED STATES PATENT 



Inventor(s): Antonio Trias Bonet 

Calle Santa Teresa 69 A 
Sant Cugat del Valles 
Barcelona SPAIN 81 190 
SPANISH CITIZEN 



Invention: SYSTEM AND METHOD FOR MONITORING AND MANAGING 
ELECTRICAL POWER TRANSMISSION AND DISTRIBUTION 
NETWORKS 



Whiteford, Taylor & Preston L.L.P. 
Seven Saint Paul Street 
Baltimore, Maryland 21202-1626 
Telephone: (410) 659-6402 



SYSTEM AND METHOD FOR MONITORING AND MANAGING 
ELECTRICAL POWER TRANSMISSION AND DISTRIBUTION NETWORKS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 This application is based upon and claims benefit of copending and co-owned U.S. 

Provisional Patent Application Serial No. 60/424,351 entitled "Method and System for 
Monitoring and Managing Electrical Power Transmission and Distribution Networks", filed with 
the U.S. Patent and Trademark Office on November 6, 2002, by the inventor herein, the 
specification of which is incorporated herein by reference. 
10 BACKGROUND OF THE INVENTION 

1. Field Of The Invention. 

The present invention relates to monitoring and management of electrical power 
transmission and distribution networks, and more particularly to a system and method for 
determining the grid state and transmission line capacity of such a network by determining the 
15 network load flow using a deterministic, non-iterative, real time analysis of the network. 

2. Description of the Background. 

The global electric industry is facing a number of challenges: an aging infrastructure, 
growing demand, and rapidly changing markets, all of which threaten to reduce the reliability of 
the electricity supply. Currently, deregulation of the electricity supply industry continues, 
20 although somewhat more cautiously than before due to California's recent experience. 

Deregulation and the drive to increase efficiencies in power systems have been particularly 
relevant in the attempt to develop new processes for intelligent observation and management of 
the grid. 
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Increasing demand due to economic and demographic variations, without additional 
generation investments, has led transmission and distribution systems worldwide to their limits 
of reliable operation. According to the North American Electric Reliability Council (NAERC), 
transmission congestion is expected to continue over the next decade. Growth in demand and the 
5 increasing number of energy transactions continue to outstrip the proposed expansion of 
transmission system. In the same line, the Edison Electric Institute indicates that the U.S. 
transmission system requires nearly $56 billion in new investment over the next decade, but only 
$35 billion is likely to be spent. Figures from the Federal Energy Regulatory Commission 
(FERC) place the total transmission congestion costs nationwide at several hundred million 
10 dollars. 

One action FERC is taking to improve coordination on the US grid is to create Regional 
Transmission Organizations (RTOs). Yet, even this important step towards nationwide 
coordination raises concerns about transmission reliability. In its report, "Reliability Assessment 
2001-2010," the NAERC stated, "The transition period from existing grid operation 

15 arrangements to the new world of RTO-managed grids may create some negative system 

reliability impacts. New system and organizational structures will need to be implemented over 
very aggressive time lines." Furthermore, the Transmission Rights market is just beginning. In 
the US FERC, as a result of three conferences, issued a working paper where the important 
characteristics were defined: LMP (Location Marginal Pricing) as the system for congestion 

20 management, the availability of a non-discriminatory standard "Network Access Service," RTO 
operation bid based day ahead and spot markets, holder's ability to sell transmission rights, and 
mitigation through market bidding rules. 
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Therefore, today more than ever before, the need exists for adequate methods for 
determining the basic functions that provide System Operators and Regional Transmission 
Organization managers with the best knowledge on their existing grid. Tools that help reduce 
the uncertainty or "fuzzy-zone" for safety operations with accurate computation of the grid state 
5 and transmission lines capacity are therefore required. 

The primary objective of operation and security management is to maximize 
infrastructure use while concurrently reducing the risk of system instability and blackouts. One 
specific type of transmission system voltage instability is the slow spreading uncontrollable 
decline in voltage known as voltage collapse. 
10 Electricity providers try to avoid power disruption to their customers. From a momentary 

interruption to a full blackout, any disturbance is costly to the provider and consumers alike. Six 
days of rolling blackouts in 2001 cost Silicon Valley businesses more than $1 billion according 
to the San Jose Mercury News. A report released by the Electric Power Research Institute's 
(EPRI) Consortium for Electrical Infrastructure to Support a Digital Society (CEIDS) notes that 
15 U.S. businesses lose over $45 billion annually from outages. 

The electrical power network is represented through the power system model by means 
of an accurate representation of all of its components: bus bars, lines, transformers, loads, 
generators, DC couplings, shunts, etc. These elements are modeled using a set of values defining 
its state (voltage, angle, and active and reactive power for nodal elements and complex flows for 
20 link elements). These values are not independent. They must satisfy the Ohm and Kirchov 
Laws, which for these variables becomes a system of non-linear equations. 

This system of equations well known as the Load Flow or Power Flow equations can or 
cannot have a solution (Voltage Collapse) and the mathematical solution to this problem 
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normally is multiple, with a degree of multiplicity as high as 2 N where N is the number of buses 
in the network. From this set of 2 N solutions, only one corresponds to the physical situation. 
The rest of the solutions are spurious and cannot represent the physical solution of a real power 
system. A standard approach to this highly nonlinear problem has been the use of numerical 
5 approximation methods. 

The topology of the actual representation can vary if the model is only detailed up to bus 
bar level, which may suit off-line studies for Planning Engineers. Yet for operations, the model 
must reach switching levels. Modeling for other purposes can also be done, as described in U.S. 
Patent No. 6,202,041 to Tse et al., which discloses a modeling method for small perturbation 
10 stability, as well as U.S. Patent No. 6,141,482 to Flint et al., which discloses an AC power line 
network simulator. 

Real time instruments in the field measure some of these parameters that are sent through 
communication lines to centralized control centers. SCADA (Supervisory and Data Acquisition) 
Systems are the basic hardware-software basis for observation and operation of a power system 

15 network (alarms, Automatic Generation Control or "AGC," etc.), and EMSs (Energy 

Management Systems) include more advanced software applications which implement the 
process of information transformation within such control centers calculating load flow, optimal 
power flow, contingency analysis, etc. For example, U.S. Patent No. 5,181,026 to Granville 
discloses a system for measuring voltage, phase angle and line temperatures in power lines. 

20 A power system model with a complete set of exact measurements for all parameters is 

not possible; hence, observation of real values is limited to a subset of all needed parameters. 
The remaining values must be estimated. Therefore, to a given set of real time values at an 
instant t are added the corresponding complementary estimated values. In order to represent a 
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feasible electrical state of the power system, these values must satisfy the Load Flow equations. 
Hence, at the heart of any real time system modeling lie two basic processes: state estimation and 
load flow equations solving methods. 

Most state estimation methods today define an external model (being the neighboring 
5 power systems' topology and values) and propagate voltage values to the internal model that of 
the given power system. It is a least square function minimization process of the differences 
between the real measured values and the estimated values. 

The standard methodology for solving the load flow equations problem has been to use 
the Fast Decoupled Newton-Raphson (FDNR) algorithm. This methodology presents two majors 
10 drawbacks: 

a) Even in the case where there is a solution, FDNR may not be able to find 
it, due to the fractal nature of the convergence region of this algorithm. This is inherent to the 
iterative nature of the Newton-Raphson Methodology. 

b) FDNR cannot assure that a solution (one that solves the mathematical 

15 equations) really represents the physical one. Newton-Raphson can jump from the neighborhood 
of one solution to the neighborhood of another in an uncontrollable way. 

The problems of the FDNR methodology are well known by the electrical sector, taking 
the form of stochastic non-convergence or dependency of the result in the order of the actions 
over the network. 

20 Several attempts to overcome these difficulties have been undertaken in the past, but with 

limited success. For example, load flow and state estimators currently used in electrical 
advanced applications at control centers, represent the state-of-the-art technology: Newton- 
Raphson Iterative methodology, as well as variants for improving convergence and speed of 
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computation (Fast decoupling, etc.), avoiding triangulation of the Jacobian, as well as new 
approaches using fuzzy logic and genetic algorithms. The list of references on this matter is not 
exhaustive but its length indicates that it is a problem yet to be solved to complete satisfaction. 
Once the model of the power system is validated as an accurate one (model topology has 
5 been improved and quality of measurements has been attained or at best ranked adequately), 
through state estimation and load flow calculation, many other processes typically take place 
within an EMS operator working environment, including: 

1) Limit violation control of parameters outside operating limits. These processes 
may comprise intelligent methods that generate proposed remedial actions by means of using 

10 load flow on the last estimated snapshot or state of the power system by the EMS State Estimator 
automatically (by means of an algorithm) or manually using a real time network simulation by 
the operator. Some physical devices, such as protections and others with or without local 
intelligence, have also been developed as alternatives, including U.S. Patent No. 5,428,494 to 
Ahuja, which presents a system for over-voltage and under-current protection, and U.S. Patent 

15 No. 5,327,355 to Chiba et al., which presents a fuzzy logic, based method for tap transformer 
settings for voltage control. Extreme remedial action always involve load shedding, which 
process is treated in some form in U.S. Patent Nos. 4,324,987 to Sullivan, II et al, 4,337,401 to 
Olson, 4,583,182 to Breddan, and 5,414,640 to Seem. A method for controlling voltage and 
reactive power fluctuations in adjacent power systems is discussed in U.S. Patent No. 6,188,205 

20 to Tanimoto et al. 

2) Planned maintenance outages assessment through instant real time on line 
simulation from the actual network state. 
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3) Optimal power plow for objective functions such as losses minimization through 
reactive power cycling. 

4) Voltage stability analysis, which can be viewed as the aggregation of the 
following: 

5 a. PV and QV curves construction. 

b. Determination of voltage collapse point and current operating point as 
well as voltage stability criterion. 

c. Generating a metric to voltage collapse. One such example is the margin 
to voltage collapse defined as the largest load change that the power system may sustain at a set 

10 of buses from a well defined operating point, as described in U.S. Patent No. 5,745,368 to Ejebe 
et al. 

d. Voltage stability assessment and contingency analysis and classification. 
Concerning voltage stability security assessment, state of the art load flow methodologies in 
general do not work. A well detailed explanation on which of these tasks they tend to fail can be 

15 found on U.S. Patent Nos. 5,594,659 and 5,610,834 to Schlueter. Because of this, Newton- 
Raphson is ill conditioned for the situation. In the cited patents, Schlueter states that current 
methods lack diagnosis procedures for determining causes of specific voltage instability 
problems, as well as intelligent preventive procedures. His method is an attempt to overcome 
this situation in certain cases. He provides for detecting if certain contingencies (line outages 

20 and loss of generation) related to reactive reserve basins can cause voltage instability. 

Another approach is that given in U.S. Patent No. 5,642,000 to Jean-Jumeau et al. where 
a performance index is related to the load demand and not to voltage for the first time. This 
index allows for determining the amount of load increase the system can tolerate before the 
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collapse, and when collapse is to be originated by a contingency, it gives a measure of its 
severity. It overcomes the computational burden of the high non-linearity of order 2 N by 
inventing a new characteristic linear equation of the exact saddle-node bifurcation point of order 
N: "Decoupled, parameter-dependent, non-linear (DPDN) dynamic systems as ones whose 
dynamics can be represented by a set of non-linear equations with a varying parameter that can 
de decoupled from the remainder of the equation". 

A method in U.S. Patent No. 4,974,140 to Iba et al. discloses discriminating voltage 
stability from the method of multiple load flow solutions. 

Also, U.S. Patent No. 5,745,368 to Ejebe et al. compares three approaches to determining 
an alternative voltage collapse point and an index, using a comparison of the method introduced: 
the Test Function Method (TFM) with two other prior art existing methods, namely, 
Continuation Power Flow (CPF) and Multiple Load Flow Method (MLF). 

Other approaches that are innovative yet still inefficient include those of U.S. Patent No. 
5,629,862 to Brandwajn et al. using artificial intelligence rule-based systems, or U.S. Patent No. 
5,625,75 1 to Brandwajn et al. for contingency ranking. 

e. Future near-term dynamic voltage stability. One such example for a near 
term definition of 25 minutes is U.S. Patent No. 5,796,628 to Chiang et al. where system voltage 
profiles are predicted and loads and contingencies are analyzed on this near-term scenario in 
terms of load margins to collapse. Continuation load flow technique CPFLOW (predictor 
corrector type of continuation power flow with a step-size control) through the nose of PV QV 
curves (saddle-node bifurcation) is reported to work without numerical difficulties. Yet, the 
patent preferred embodiment describes the sensitivity of the number of final iterations to the 
attainment of a good approximation point for the next solution by the predictor. It is also stated 
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that good step-size controls are usually custom-designed for specific applications. So again, 
there is some craftsmanship as in all PV QV curve construction using any derived method from 
Newton-Raphson iterative process. 

5) On-Line transient stability. This is a more ambitious task, entering the realm of 
5 the differential equations where the right hand term is a load flow equation. U.S. Patent No. 

5,638,297 to Mansour et al. defines via an artificial contingency on-line transient stability 
assessment. 

6) Load forecast. We list here this process even though it is not related to load flow 
methodologies because knowing the forecasted load profile will help in many instances while 

10 analyzing future contingencies and generating action plans (limits back to normal, restoration). 
Standard methodologies used by successful methods include the more classical autoregressive 
methods (ARIMA) Box Jenkins time series approach, as well as more recent artificial 
intelligence neural network approaches. 

7) Disturbance detection and restoration. 

15 a. For distribution grids, the problem is more simple and well known. 

Restoration can be managed through a set of rules (small expert systems) since the topology is 
radial. State of the art is mostly centered in fault location and its resolution as well as protection 
schemes by different standard and creative ways. Patents include U.S. Patent Nos. 5,303,1 12 to 
Zulaski et al., 5,455,776 to Novosel, 5,568,399 to Sumic, 5,684,710 to Ehlers et al., 5,973,899 to 

20 Williams et al., 6,185,482 to Egolf et al., and 6,347,027 to Nelson et al. 

b. For transmission grids, the restoration problem has not been solved 
satisfactorily as a general universal procedure valid for any power system network. With ageing 
infrastructures and growing demand, disturbances are increasingly likely to happen. 
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Traditionally, restoration after a disturbance has been one of the most difficult things for 
electrical companies to handle. While hundreds of hours of systems analysis and documentation 
go into restoration plans, they never match the reality of any specific disturbance, and they are 
dynamic in nature. Automatic tools for helping operators have been attempted. Avoiding the 
need for local rules specification would be desired. Detection is related to intelligent alarm and 
topology changes processing. Restoration plan validation by the operator requires load flow 
calculation after each step in order to guarantee a feasible electrical network state after each and 
every action, with the post-disturbance steady state as initial condition of the action plan. 

As we have seen, all of the above central processes need a working, real-time load flow 
method. These methodologies on which the industry has based, up until now, the on-line real 
time monitoring and managing of networks as well as off-line analysis tools for planning, 
programming, and for investing decisions support, generally behave well under certain 
continuity of the network condition. Iterative in nature, they need initial points near the solution 
or equivalent knowledge of the previous solution to compute the next stage in a real time 
environment. This last aspect is responsible for not being able to behave well when disruptions 
of the system state take place, when a major disturbance or blackout takes place. To conclude, 
we add that when the electrical network state is close to voltage collapse, precisely when 
operators and planners need the support of these tools the most, traditional methods fail and 
frequently cannot deliver a correct calculation. 

SUMMARY OF THE INVENTION 

Disclosed herein is a deterministic non-iterative method that improves the existing 
methods to solve the load flow equations of any power system. Such method in turn provides 
improved methods for state estimation, generation of restoration plans, the construction of PV 
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and QV curves, voltage stability and contingency analysis, optimal power flow, and operation 
limits control. 

In a preferred embodiment of the method of the invention, a physical solution of the 
central load flow problem is found using the following steps: 

a) Embed the load flow problem in a parametric homotopy that goes continuously 
from the 0-case to the problem case; 

b) Develop in power series the values of the equation's unknowns in the 
parameters) of the homotopy in a neighborhood of the 0-case value of the parameter; and 

c) Use analytical continuity to find the value for the equation's unknowns in the 
problem case. 

For suitable analytical continuation techniques using algebraic approximants (continued 
fractions, for instance), the above-described procedure always gives the correct solution (i.e., the 
physical one) when it exists. If no solution exists, then we are at the voltage collapse state of the 
power system. The present invention thus relates to a constructive method for finding such a 
solution (or determining that no solution exists and thus that the system is at voltage collapse), 
and a system for employing such method. 

DESCRIPTION OF THE DRAWINGS 

Figures la-le are schematic representations of convergence regions realized through 
implementation of a prior art FDNR method on a two-bus network. 

Figure 2 is a schematic representation of a method for determining power series 
coefficients for voltages V[n]. 
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Figure 3 is a schematic representation of a method for evaluating an n-order approximant 
of a continued fraction approximation for the power series coefficients produced by the method 
depicted in Figure 2 to provide a solution to the load flow equations (L). 

Figure 4 is a schematic representation of a method employing a loadflow determination 
method of the invention for purposes of performing state estimation. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

To illustrate the well known fact of the ill condition of the existing state of the art 
methodology in the vicinity of voltage collapse, we will use a very simple example of a very 
small network with final quadratic load flow equations (exact solutions are easily computed). 

The general load flow problem has at least one swing node and a set of nodes (generation 
and/or load). In the very simple minimalist example chosen, we will only have one node as well 
as the swing node. This swing node does not vary its voltage value no matter how load and 
topology may vary in the rest of the network. It acts as a large generator or substation capable of 
providing any power required by the system. Only active and reactive power are calculated at 
the swing bus, balancing the sum of both at the rest of the nodes. 

An alternate current in stationary regime satisfies Ohms law with complex values. This 
is the origin of the complex values (X inductance, V voltage, R reactance) used in our problem. 

Ohms Law for this simple system is simply: 

V-VO = Z I 

where VO is the initial voltage, I is the complex intensity, and Z the impedance. The trivial 
solution is 

V= VO + ZI 
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Since VO and Z are known, if the intensity I consumed at the load node is given, the complex V 
value is easily computed (singular Z=0 is short circuit and we exclude it). 

In our example, the circuit has only 2 nodes or buses: the swing one with complex 
voltage fixed at VO = 1, i.e., |V0| is 1, and its angle or phase is 0. 
5 The other node is a load, and the goal is to calculate the value of the complex voltage: 

module and phase equivalent to real and imaginary parts. 

The network has only one link joining both nodes with an impedance: 

Z = R + j X (j complex unit, R Reactance, X inductance). 
In general you do not know I (not easily measurable values in high and medium voltage nodes), 
10 which could reduce the load flow problem to a linear one easily solved by matrix inversion. 

Loads are only known at complex power values, that is P: Active power and Q: Reactive power, 
the first value being the object to be billed, and easily measured at transformers level. Let S be 
the complex power: 

S = P+jQ 
15 The relation among these is: 

S = VI* 

Where * stands for the complex conjugate. Therefore, Ohms Law becomes: 

V = V0 + ZS*/ V* 

which unfortunately is the quadratic and non-linear equation that has to be solved for larger N: 
20 number of nodes. This is the real difficult problem of Load Flow calculation. 

In this simple example, the equation system may be solved as follows. Solution for 
VOol is similar to the case V0=1 . Let H be 

H = Z S* = (RP + XQ) + j (XP-RQ) 

V V* = V* +H 



VR2 + Vl2 = VR - j VI +HR +j HI 
The imaginary equation gives us VI: 
VI = HI 

5 

Substitution in the real part: 

VR2 + ffl2 = VR + HR 
VR2-VR-(HR-Hl2) = 0 
VR = V 2 +- sqrt(l/4 + HR - ffl2) 

10 

We also obtain the Power real part. 
HR - Hl2 < -1/4 there is no physical solution. 

HR - Hl2 > -1/4 the solution is double: one physical and one spurious (i.e., spurious is equal 
to Vf + 2k n where Vf is the correct physical solution). 

15 

HR - Hl2 = -1/4 in the limit both solutions coalesce, both branches coincide and with more load 

there is voltage collapse. 

Figures la through le provide schematic representations of convergence regions realized 

through implementation of a prior art FDNR method on a two-bus network. More particularly, 
20 the figures show a two bus network with one swing bus. Parameters of the basic example 

depicted in Figure 1 are: R=0, X=0.2, P=-0.8 and Q=-0.2. In Figures la through le, the Reactive 

Q will be varying down to -1.15 approaching voltage collapse. 

The legend for Figures la through le is the following: initial points for the voltage 

(module and angular representation) that lead the load flow state-of-the-art Newton-Raphson 
25 Fast Decoupling methodology to the physical solution are colored white, while the black color 

depicts those points where the Newton-Raphson Fast Decoupling methodology would produce a 

spurious solution, a non-physical solution, and when there is no convergence (and thus no 

solution) even though we know that a solution exists. 
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Figure lb shows how the simple initial values of the model, within a region (white) of 
initial points, allow for the iterative process of the state-of-the-art methodologies to converge to 
the solution. Still the voltage collapse is not near, since the spurious solution S is far away from 
the physical one P and most of the remaining area comprises points that if used as initial ones 
5 will allow the Newton-Raphson Fast Decoupling method to converge to this spurious solution. 

Yet in Figure lc, we slowly approach reactive power parameters that makes the system 
reach the voltage collapse zone, and each time we continue this procedure we see a behavior of 
the state-of-the-art methodologies that worsens when nearing voltage collapse. 
As shown in Figure Id, in which the black and white areas are interspersed, initial points colored 
10 black that are used as starting points will produce incorrect solutions being in the boundary of the 
white regions. This can lead to misleading results. More specifically, the situation is even worse 
because the spurious solution area is intermixed with the convergence area. 

Finally, in Figure le, we observe that in spite of the existence of solutions, the state-of- 
the-art methodology can only find them at very few initial points within the black area since at 
15 this value for the reactive power Q, the system is very close to the voltage collapse attained at 0.5 
as the solution to the load flow equations. 

In this simplistic example, the new method to be introduced behaves in an excellent 
manner with regard to approaching voltage collapse. For this problem, all the area would always 
be colored white (indicating that every point would lead to the physical solution). 
20 Extending the chaotic behavior of this two bus model to larger networks, it is clear that 

unreliable results can be introduced near voltage collapse for transmission grids. 

The method of the invention is a deterministic, non-iterative process to finding the 
solution to the load flow problem that behaves well near voltage collapse. The method 



converges universally if the problem admits a solution, and never if the problem does not have a 
feasible physical solution. The following discussion provides a constructive procedure for 
producing such solution to the load flow problem. However, in order to present such procedure, 
it is first necessary to establish the following principles: 
5 1 . The physical solution must be connected in a continuous way to the non-load and non- 

generation case (0-case), in which all the voltages are equal to the nominal voltage, and there is 
no energy flow in the links. The reason for this lies in the fact that the 0-case is physical (it is 
possible to build a real electrical power system with this state) and any other physical state can 
be reached by increasing simultaneously in a continuous way, load and generation until the final 

10 state is reached. 

2. The quantities that appear in the equations (voltages, power, and flows that are 
complex numbers) are constrained to have functional relations between them with a very strong 
property called analyticity. This is a property of functions defined in the complex plane that 
reflects deeper symmetries of the system than is represented by the functions. In this case, 

15 analyticity is a property implied in the definition of the Ohm and Kirchov laws and, thus, by the 
load flow equations. 

Using these two facts as a framework, we define the method for finding a physical 
solution of the load flow problem in the following steps: 

a. Embed the load flow problem L in a parametric homotopy L(s) that goes continuously 
20 from the 0-case (L(0)) to the problem or objective case (L(l)). 

b. Develop in power series the values of the equation's unknowns in the parameters of 
the homotopy in a neighborhood of the 0-case value of the parameter. 



c. Use analytical continuity to find the value for the equation's unknowns in the problem 

case. 

For suitable analytical continuation techniques using algebraic approximants (continued 
fractions, for instance), it is possible to prove that this procedure always gives the correct 
5 solution (i.e., the physical one) when it exists. If no solution exists then the power system is 
undergoing voltage collapse. 

Details of the basic steps to calculate the solution with the method of the invention for 
general N are the following. 

First, we construct the embedding transforming the load flow equations into a function of 
10 a single complex variable. 



-> Us) 



L(s) 



S" S* 

15 + y^ 2 + - + y,v v n = -r y^Ms) + y l2 v 2 (s) + ... + y lN v n (s) = s-± 



. _ s: 



y« v i + y 2 2 v 2 + - + yis^n = -4 y 2 Ms) + y 22 v 2 (s) + ... + y 2N v„(s) = s-f 



V 2 V 2 



s* s' 

j>m v > +y N 2^2 +-+y NN v« =-r- y N Ms)+y N 2v 2 (s)+...+y NN v n (s) = s^- 

20 



For an n-bus case, let Yij be the admittance matrix of an n-buses network (0 is a swing 
bus), with Si and Vi the complex power and complex voltage at bus i. The loadflow equations 
can be written as 

25 

|r = 2%F* (0.1) 

"i *=o 
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In order to solve the load flow equation, we define an embedding in a family of problems 
depending on a parameter s such that we know the solution for s=0, and for s=l we recover the 
5 original equations. One of the possible embeddings is: 



10 



20 



v i \ s ) k=0 k=0 



V 0 (s) = l + (l-s)V 0 

v k (0)=iyk 

V k (l) = V k ,\/k 



(0.2) 



Next, we define a functional transform from the analytical functions to the infinite sequences set: 



r:/00->/[n] = l 



d"f(s) 
ds" 



s=0 



where f[n] is the n coefficient of the MacLaurin series expansion of f(s). 

A*) = /[0] + /[l> + /[2]5 2 +...+f[nV+... 

15 with the properties 

r(f(s)) = f[n] 
<sf(s)) = f[n-l] 



r(f(s)g(s)) = (f*g)[n] = ^f[k]g[n -k] 



k=0 



We rewrite (0.2) as 



k=\ 



= sSXis) - Y i0 V 0 (s) - (1 Y ik 

1 



And applying the transform to both sides of the equation we get 



(0.3) 



(0.4) 



(0.5) 



(0.6) 
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&/ 1 w=^ , [»-i]-^(i-ft.-^Fo)-(^-^)a l 



(0.7) 




defining a recurrence over n taking into account that 



Wi{s)Vi(s) = 1 => {W t * V,)[n] = 6, 




(0.8) 



W i [n] = -£ t W i [k]V i [n-k] 



* being the sequence convolutions operator. 



The steps to calculate the coefficients in the series expansion to order n, are 



i) 



Initialization 



V i [0] = W i [0]^\ 



(0.9) 



ii) For m=l to n 

Calculate Vi[m] solving the linear system (1.7) 
Calculate Wi[m] using (L8) 
The entire process is represented in Figure 2. 

This will give the power series expansion of Vi(s) up to order n. In general, however, 
this series will not converge for s=l. Nevertheless, a continued fraction expansion of the power 
series will converge for all s values when voltages are given inside the solution set continuously 
connected to the s=0 case (no load). 

Next, from the series coefficients, it is possible to build a rational approximant for the 
function obtained by analytic continuation from the point s=0 to s=l. There is a proof assuring 
that if the set of equations has a solution in the physical branch, it is always possible to find a 
continuation path from s=0 (no charge) to s=l, free of singularities, and obtain the solution to the 
equation by evaluating the rational function for s=l. 
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An algorithm for constructing an algebraic approximant (e.g., a continued fraction 
approximation) to a power series is the well known Viskovatov method, as described in A. 
Bultheel, "Division Algorithms for Continued Fractions and the Pade Table," J. Comp. Appld. 
Math. No. 6, 1980, which is incorporated herein by reference. Another methodology is to use 
Pade-Hermite Approximants or any technique capable of computing an algebraic approximant 
from the power series of an analytical function, as described in George A. Baker and Peter 
Graves-Morris, "Pade Approximants, Second Edition," Encyclopedia of Mathematics and its 
Applications, Volume 59 (Cambridge University Press, 1996), which is incorporated herein by 
reference. 

For clarity, we will explain the Viskovatov approach used within our process. It is 
inspired in a "double inversion" of the power series. 
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m = no] + m + f[2]s 2 + ... + f[nV + ... = 

= /[O] + 5(/[l] + /[2]» + ... + /[«>"- +...) = 

5 



= /[0] + - 



1 



= /[0] + 



/[l] + /[2>+...+/[n> n -'+. 

5 



/ (,) (*) 

with 

f (l \s) = r- 

f[l] + f[2]s+...+f[n]s'- l +. 

/ (1) [0] + / (1) [1> + ...+ f 0) [n -\]s"-* + .. 



f(s) = f[0) + - 



/>] + - 



1 



= /[0] + - 



/ (1) [l] + ...+/ (,) [»-l>"- 2 +... 

s (0.10) 



/ 0) [0] + 



/ (2) [0] + - s 



/ (3) [0]+... 

The f (,+1) power series calculation can be performed using the i* 0 power series applying 
the (0.9) recursion set forth above. Here, the power series f particularly corresponds to the v 
function (voltage) and w (its inverse). 

After that, the n-order approximant An(s)/Bn(s) of the previous continued fraction can be 
evaluated using the recursion relation 

Ms) = / (i) [0]4-.(s) + sA i _ 2 (s),i = 2,3,4... 



B 0 (s) = l,B l (s) = f (l) [0) 

W = /®[0P H (j) + rff M (*),' = 2,3,4.... lU ' 
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Evaluating the n-order approximant An(s)/Bn(s) for the calculated Vi(s) power series in (0.9) for 
s=l, giving V(s=l) will give the solution to the original loadflow problem, as can be seen from 
(0.2). 

Figure 3 shows the scheme of the computational process to get the n-order approximant 
for the calculated Vi(s). More particularly, we begin from the power series coefficients for the 
voltages V[n] (calculated using the schema of Figure 2). Using those power series coefficients, 
and applying an algebraic approximant (e.g., a continued fraction methodology, such as 
Viskovatov, Pade-Hermite Approximants, or other continued fraction methodologies), we build 
the f[n] continuous fraction. In order to evaluate, it is necessary to build the series of 
approximants A[n]/B[n] using the f[n] and the two previous coefficients of the A and B series. 
Finally, using the approximant of high order A[n]/B[n], and evaluating them for s=l, we are, in 
fact, calculating the V solution to the initial problem. 

The above described method for determining the loadflow in an electric power generating 
system may be employed in a number of aspects for general management of the electrical grid, 
including observation and estimation of the network state, voltage stability and contingency 
analysis, optimal power flow, limit controls, and system restoration following a voltage collapse. 
Described below in greater detail is one such aspect concerning state estimation. However, other 
applications may likewise utilize the above-described method, particularly including: the 
generation of dynamic restoration plans as a path search method; generalized OPF as a path 
search method and limit controls as a boundary case; improved methods for generating PV and 
QV curves indirectly through substitution of available load flow techniques by the above- 
described method; determination of voltage collapse region characteristics using zeroes and 



22 



poles of the approximants; and voltage stability analysis and contingency analysis indirectly 
through substitution of available load flow techniques by the above-described method. 

With particular reference to a method for state estimation, reference is made to Figure 4, 
which provides a schematic overview of such method. Data coming from the field includes 
loads generations, voltages, flows, and the state of breakers, among others. The network 
topology describes the possible connectivity of the electrical network. 

After receiving the data from the field and the topology of the model, a battery of tests 
1 102 using Artificial Intelligence are done in order to make inferences on the missing 
information and quality of the available measurements from the electrical laws. The tests 
include logical considerations about the connections and the measures observed. 

Dynamic assignment 1 104 of a quality parameter which is historically followed is done, 
such quality parameter being the result a very robust estimator with no need for a very high 
percentage of observable measures. The quality parameter is calculated from the historical 
comparison between the field measurement and their estimated value. This quality parameter 
expresses the confidence in the field measurements and is used as a weight in the estimation 
process. It calculates for the complete network, thus avoiding the effort traditionally necessary 
to work with an external model and then propagate to an internal one. 

The State Estimation process consists on standard least square minimization on the 
weighted differences and takes place using Gauss Seidel. 

The above-described method for determining the loadflow in the network allows us to 
accept only feasible physical states (continuously connected to the no-charge solution). Only 
these states can be seen in the field. Hence, we require every state estimation to always run the 
loadflow 1108. 
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If we do not get a solution and the electric system is not at voltage collapse 1110, there 
can be no more than three reasons: (i) synchronization problems of snapshots 1112 (measures 
from different temporal intervals are used); (ii) measurement problems (error in the measurement 
device or the communication line); or (iii) modeling problems (errors in the static parameters 
5 defining the model) 1 1 14, 1 1 16. This simple result has allowed a very powerful diagnostic and 
calibration kit for quickly attaining an improved model closer to the physical grid than has been 
available from prior art systems, as well as improvements in the detection and, hence, correction 
of the bad quality of certain measures. We can estimate even in the region of voltage collapse 
giving utility operators the support of reliable calculations when such information is needed the 
10 most. 
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